Induced interactions and the superfluid transition temperature in a three-component 

Fermi gas 



O 

(N 

(N 



J. -P. Martikainen,^ J. J. Kinnunen,-^ P. Torma,^ and C. J. Pethick^'^ 

'NORDITA, RoslagstuUsbacken 21, 106 91 Stockholm, Sweden 
^Department of Applied Physics, Helsinki University of Technology, P.O. Box 5100, 02015 HUT, Finland 
'^The Niels Bohr International Academy, The Niels Bohr Institute, 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
(Dated: January 2, 2010) 

We study many-body contributions to the effective interaction between fermions in a three- 
component Fermi mixture. We find that effective interactions induced by the third component 
can lead to a phase diagram different from that predicted if interactions with the third component 
are neglected. As a result, in a confining potential a superfluid shell structure can arise even for 
equal populations of the components. We also find a critical temperature for the BCS transition 
in a ^Li mixture which can deviate strongly from the one in a weakly interacting two-component 
system. 
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By using Feshbach resonances to change the effec- 
tive interaction between ultracold atoms several groups 
have probed the crossover from the Bardeen-Cooper- 
Schrieffer (BCS) superfluid to a Bose-Einstein conden- 
sate of molecules |ll-[lO|. Studies of the crossover have 
provided important insights into fermionic superfluids 
around the unitarity limit of strong interactions. Im- 
portantly, it has become experimentally feasible to study 
also more complicated mixtures than Fermi gases with 
two different atomic internal states. Bose- Fermi mix- 
tures ll|, and Bose-Einstein condensates with many 



components have been created using many different se- 



tups 



14 1 . Also, heteronuclear Fermi- Fermi mix- 



tures [15|, 16 1 and even heteronuclear Fermi-Fermi-Bose 
mixtures [l7| have been recently demonstrated. In yet 



another breakthrough a three-component Fermi mix- 
ture of atoms in the three lowest hyperfine states of 
^Li 18, UM has also been demonstrated. Such multi- 



component systems have some intriguing similarities with 
quark matter counterparts where color superconductivity 
may appear 

The purpose of this Letter is to explore how the in- 
duced interactions due to the third component mod- 
ify the expected behavior of three-component mixtures. 
This is important because, depending on parameters, the 
many-body effects can change the effective interaction 
between atoms substantially and on occasion even change 
the relative magnitudes of couplings between different 
components. Such changes imply that phase diagrams 
predicted using only two-body scattering properties can 
be incorrect. Also, even when the corrections to the effec- 
tive interactions are weak, they can cause large changes 
to the critical temperature for the BCS transition. In- 
deed, for a two component system Gorkov and Melik- 
Barkhudarov (CM) showed [2l| that the perturbative 
correction to the effective interaction can reduce the crit- 
ical temperature by a constant factor of (Ae)^^^ ~ 2.22 
in the weak-coupling limit. Also in spin-density imbal- 
anced systems such corrections have been shown to have 



a considerable effect [221 . Here we will analyze the effects 
of analogous corrections on the three component system 
and find important changes to the CM result. As an in- 
teresting consequence of these many-body corrections we 
predict in a spatially varying confining potential (typi- 
cally harmonic trap) the appearance of superfluid shell 
structures even in the absence of population imbalance 
(polarization) of the components. These shell structures 
are due to many-body effects only and therefore funda- 
mentally different from earlier predictions of shell struc- 
tures due to population, mass, or trapping potential im- 
balance 2^, 24 1 We also point out that many body effects 
due to the third component provide a new way to tune 
the effective interaction between the two other fermions 
and that this contribution can dominate over the usual 
CM contribution. 

Earlier, intriquing results have been found experimen- 
tally for the critical temperature of iron-based multiband 
superconductors 25] and degenerate three-component 
Fermi gases have been studied theoreticall y i n a lat- 
Furthermore, pai ring (23l. [28l. |29| , stabil- 
and breached pairing [3l| have recently been 
studied in a three-component fermionic mixtures. How- 
ever, these theoretical approaches did not consider situ- 
ations directly relevant to ongoing experiments and also 
did not study how the many-body effects due to the pres- 
ence of the third component influence the properties of 
the other two components. Some aspects of the many 
flavor problem were discussed by Heiselberg et al. [S^ . 

The relevant second order diagrams which give rise to 
induced interactions between fermions of type 1 and 2 
are shown in Fig. [T] In these diagrams the arrows are 
the component propagators and dashed line is a contact 
interaction with strength Ua^ between components la- 
beled by a, /3 g {1, 2, 3} {a ^ (3). These couplings can be 
expressed in terms of the scattering lengths aap through 
Uap — 2'Kh^aai3/rna/3, where ruap = {l/rria + l/mp)~^ is 
the reduced mass. Of the diagrams shown, the diagram 
(a) is relevant in the case of a two component system 
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FIG. 1: Diagrams of second order in the interactions for the 
induced interaction between components 1 and 2. Sohd lines 
represent atoms and the dashed lines interactions between 
them, and numbers indicate the fermionic component. The 
interactions are taken to be antisymmetrized with respect to 
interchange of spins and consequently interactions between 
atoms of the same species are absent. 



with a contact interaction between unlike fermions [33| 
and for equal mass fermions gives rise to the GM cor- 
rection mentioned earlier. In a three-component system 
the diagram (b) describes the induced effect of the third 
component. Similar loop diagrams with the mediating 
fermion in component 1 or 2 are forbidden in the s-wave 
scattering channel for symmetry reasons. More formally 
the diagram (a) indicates the induced interaction 



/Ki(k + q/2)]-/fc(k-q/2)] 
a(k + q/2)-e2(k-q/2) 



where q 
tion 



V 



3c 



p + p' and (b) describes the induced interac- 
/[6(k + q72)]-/[6(k-q72)] 



(p,p') = C/l3C/23^ 



C3(k + q'/2)-C3(k-q'/2) 



with q' = p — p'. In these formulas ^a(k) = Ti^k"^ /2ma — 
fia are the free atom dispersion relations and /(e) is the 
Fermi distribution. 

In the weak coupling limit the scattering processes 
around the Fermi surfaces dominate and to find the effec- 
tive coupling the induced interactions are averaged over 
the Fermi surfaces. In this way we find that the effective 
coupling between fermions of types 1 and 2 becomes 
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where kp^a is the Fermi wavevector for the component a 
and we have assumed that fermions 1 and 2 both have a 
mass mi while the third component has a mass . The 
function F{y) is given by the integral 
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FIG. 2: Scattering lengths in units of the Bohr radius between 
^Li atoms as a function of magnetic field. (The figure is taken 
from Ref. fsl.) 



whose analytical solution is given by 
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The effective interactions in other channels can be found 
in the same way. 

In Eq. ([T|) the first term describes the two-body scat- 
tering in the absence of Fermi seas, the second term gives 
rise to the GM correction, and the third term describes 
the effect of the interactions with the third component 
and its Fermi sea. The correction due to the second term 
always suppresses the critical temperature for the BCS 
transition. However, the last term is proportional to the 
product a23ai3 and can have either sign. Therefore, the 
presence of the third component can either suppress or 
enhance the critical temperature. 

Since three component systems have been demon- 
strated using ^Li atoms, let us now investigate these 
many-body effects using the coupled channel scattering 
data for ^Li [3^]. In Fig. [5] we show the scattering lengths 
between different components of the ^Li mixture (|1), |2), 
|3) refer to the states iF^mp) = |l/2, 1/2), |l/2,-l/2), 
3/2,-3/2), respectively). It can be seen that, in the 
absence of many-body corrections, the 1 — 3 channel 
has the most negative scattering length for weaker mag- 
netic fields, while at magnetic fields above the Feshbach 
resonances the 1 — 2 channel eventually becomes domi- 
nant. In the simple mean- field picture one would infer 
that these channels are also the ones with highest crit- 
ical temperatures. However, when we include induced 
interactions, density dependencies appear in the effec- 
tive coupling strengths and change the simple picture in 
which the third component is neglected. 

Let us first explore the case where all ^Li components 
have the same density. In Fig. [3] we show the dominant 
coupling channel in the magnetic field-density plane. It 
can be seen that below the Feshbach resonance the 1 — 3 
channel dominates for smaller densities, but for densities 
higher than about 2 • 10^^/cm'^ there is a possibility that 
the 1 — 2 channel becomes dominant. At higher magnetic 
fields we find a possibility of dominant 1 — 3 coupling in 
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FIG. 3: The dominant interaction channel including induced 
interactions for an equal-density ®Li mixture. The solid black 
lines indicate phase boundaries. Below the (red) dashed line 
both I/cfOisI and \kpa23ai2\ are less than one, while below 
the (red) dot-dashed line both lfci?ai2l and |fc^a23ai3| are less 
than one. In (a) we show the interesting regions at low mag- 
netic fields and in (b) the behavior at higher magnetic fields. 



the region where the scattering lengths would predict the 
1 — 2 channel. In experiments the atoms are trapped, 
and applying a local density approximation (which has 
been sufficient to describe many experiments) with our 
results suggests the interesting possibility of different su- 
perfluid phases appearing in different parts of the cloud 
even in the absence of polarization or unequal trapping 
potentials/masses [23j. This possibility is a many-body 
effect caused by the induced interactions only; for a bal- 
anced system at zero temperature, simple mean field the- 
ory would not predict the shell structures that arise from 
the density dependence of the GM correction as shown 
here. 

It is important to investigate what these results im- 
ply for the critical temperature. In attractive di- 
lute Fermi gases, the critical temperature for the BCS 
transition in the weak coupling limit is ksTc/eF 
exp[-7r/(2fcF|a<^*^|)] where a<^* = U''^N{€F)TT/{2kF) is 
the effective scattering length and N{eF) is the density 
of states at the Fermi level. We now use this result 
to estimate the many-body correction to Tc in a three- 
component ^Li system. For simplicity we use the above 
functional dependence in all regions where the coupling 
is attractive, but indicate the regions where |fci?a| > 1 
in the figures. In those regions the weak-coupling for- 
mula is only suggestive. In Fig. 2] we show the fraction 
Tc/Tcfi for the equal density ^Li mixture {Tc^ is the crit- 
ical temperature in the absence of induced interactions) . 
As can be seen, the correction to Tc is often very dif- 
ferent from the 1/2.22 « 0.45 GM result and shows a 
non-trivial behavior as a function of the magnetic field 
and density due to complicated variation of the scatter- 
ing lengths. This also makes it possible that the pairing 
channel is changed due to many-body effects. At high 
fields above the Feshbach resonances it is possible that 
the critical temperature is enhanced since the effective 
scattering length there becomes more negative due to in- 
duced interactions. However, this happens in the region 
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FIG. 4: The highest fraction Tc/T^.o of the critical tempera- 
tures with and without induced effects for the equal-density 
''Li mixture. The fraction is only computed in the regions 
where the dominant effective interaction is attractive and set 
to zero elsewhere. Below the yellow line |fcFa''''^| < 1. Regions 
where the optimal pairing channel is changed are visible as 
kinks in the fraction Tc/T^fi- The inset shows a close-up into 
the region of high magnetic fields and low densities. 



of stronger interactions where our results are not neces- 
sarily quantitatively accurate. 

In Fig. Owe demonstrate another possibility for chang- 
ing the critical temperature: the use of density imbal- 
ance. In Fig. [5] (a) we show an example of how Tc in the 
1 — 3 channel is changed as the density of the compo- 
nent 2 is varied. It is again clear that the result deviates 
substantially from the GM prediction, but the Tc is nev- 
ertheless suppressed by the component not involved in 
pairing. In Fig. [5] (b) we show the similar result in the 
1 — 2 channel which dominates at higher magnetic fields. 
Due to the different behavior of the scattering lengths, 
here the induced interactions can act to enhance Tc above 
the value predicted by the usual mean-field theory. 

Finally, since heteronuclear fermionic mixtures are 
experimentally feasible 15|, |l7[ let us briefly discuss 
what our results imply in that case. A mass imbal- 



ance can be realized if the third component is a dif- 
ferent isotope, but also if the third component experi- 
ences an optical lattice which changes its effective mass. 
In the latter case, for the formulas derived here to be 
valid, the filling fraction of all the components should 
be much less than one. For higher filling fractions the 
Fermi surface is no longer spherical and the result would 
change considerably [35|]. We focus on a scenario with 
equal masses for atoms of type 1 and 2, since it is 
known that unequal mass of the interacting fcrmions sup- 
presses the critical temperature 0, and for this rea- 
son the equal mass superfluidity appears more generic. 
The effective interaction between 1 and 2 is given by 
Eq. (H)). Note that since {m^ + mi)^/4miTO3 > 1, in- 
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FIG. 5: The fraction Tc/Tc.o of the critical temperatures with 
and without induced effects for the ®Li mixture. In (a) we 
show an example of how the critical temperature is changed 
in the 1 — 3 channel when the density of the component 2 is 
varied while in (b) we show the same for the 1 — 2 channel 
when the density of the component 3 is varied. We only focus 
on those magnetic fields where the pairing channel in question 
is the dominant one. 



duced interactions become relatively stronger in an un- 
equal mass mixture and mass imbalance can be used 
to enhance the role of many-body corrections. If 6 = 
(a^^/laigaasl) {kF.i/kp^zf F{1) / FikF^/kp^s) and b > 1, 
the contribution of the third component to the induced 
interaction becomes larger than the GM contribution 
when TOa/mi > (26 — 1) + 2-\/6(6 — 1) or when m^/nii < 
(26-l)-2y&(6^nj. If 6 < 1 the contribution due to the 
third component is dominant for all mass ratios. How- 
ever, for mass ratios larger than about 13.6 other physics 
can come into play since then weakly bound diatomic 
molecules might become coUisionally unstable [36| . 

In this Letter we have explored the induced inter- 
actions and their role in the BCS pairing in a three- 
component Fermi mixture. We found striking differ- 
ences from physics ignoring these many-body corrections. 
In particular, we found that when the induced inter- 
actions are taken into account, the phase-diagram can 
change drastically, that shell structures in traps can ap- 
pear even without number, mass, or trap imbalance, and 
that the critical temperature for the BCS transition is 
strongly dependent on the induced interactions in the 
three-component systems. 

We thank Academy of Finland (Projects No. 213362, 
217041, 217043, and 210953). 
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